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Abstract 

An effective low energy field theory is developed for a system of two chains. The main novelty of 
the approach is that it allows to treat generic intrachain repulsive interactions of arbitrary strength. 
The chains are coupled by a direct tunneling and four-fermion interactions. At low energies the 
individual chains are described as Luttinger liquids with an arbitrary ratio of spin v s and charge v c 
velocities. A judicious choice of the basis for the decoupled chains greatly simplifies the description 
and allows one to separate high and low energy degrees of freedom. In a direct analogy to the 
bulk cuprates the resulting effective field theory distinguishes between three qualitatively different 
regimes: (i) small doping (v c « v s ), (ii) optimal doping (v s « v c ) and (iii) large doping (v s « v c ). 
I discuss the excitation spectrum and derive expressions for the electron spectral function which 
turns out to be highly incoherent. The degree of incoherence increases when one considers an array 
of ladders (stripe phase). 

PACS numbers: PACS numbers: 71.lQ.Pni, 72.80. Sk 



I. INTRODUCTION 

In their seminal paper Dagotto and Rice [l| came with an idea that one can get an insight 
into physics of copper oxide materials by studying just two coupled CuO chains. This gave 
rise to the belief that the two-chain problem contains in itself seeds of the rich physics of 
the cuprates. Since then the interest to the problem of ladders (that is how such two-chain 
systems are usually called) has never faded and there have been numerous attempts to tackle 
it using both numerical and analytical approaches. There are also experimental examples of 
fermionic ladder materials, such as "telephone number" compound Sr 14 _ a; Ca 3; Cu 2 404i, which 
display similarities and differences with the uniform cuprates . There is a possibility 
that a model of weakly coupled ladders may describe the cuprates with stripe ordering (for 
a review, see jj]). 

As far as the theory of fermionic ladders is concerned, the efforts to tackle it can be 
separated into numerical and analytic ones. Numerical techniques allow to deal directly 
with the experimentally relevant range of the model parameters. A fermionic ladder is 
usually modeled by the t — J model with equal tunneling matrix elements along the rungs 
and legs of the ladder, the ratio of the exchange to the tunneling matrix element is taken as 



^he-art direct 
8| respectively 



J/t = 1/3 — 1/2, the ring exchange is also considered. Examples of the state-o: 
diagonalization results and DMRG calculations can be found in 0,3 and Q, 
. The use of the t — J model assumes that the on site Coulomb interaction on the individual 
chain is significantly greater than the tunneling matrix elements. At such circumstances 
one may expect a significant difference between the spin and the charge velocities for an 
individual chain and the numerical approaches obviously take this feature into account. 



The analytical approaches follow t 
renormalization group theory (see 



re seminal papers [9] and are based on bosonization and 
and references therein). They all start from the limit 



of weak interactions and have a great problem with tackling the difference in velocities. For 
approximately equal spin and charge velocities there is a remarkable result first obtained 
in {9]. It turns out that for a quite generic choice of bare coupling constants the theory 
scales to a strong coupling state with an enlarged U(l)xO(6) symmetry greater than the 
symmetry of the lattice Hamiltonian. This feature allows one to solve the low energy model 



exactly 



llj and to calculate its correlation functions 11], 12]. The phase diagram of the 



strongly correlated model contains a phase with a d-wave superconducting quasi long range 
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order; this has been taken as evidence that such model may be relevant for the physics of 
cuprates. 

Both numerical and analytical approaches have their drawbacks. The present day numer- 
ical techniques used to study 2-leg ladders are exact diagonalization (Lanczos method) and 
DMRG. Although Lanczos method has generated many impressive results for the problem 
of doped ladders (I discuss them further in the text), it does not allow to study long chains 
nor (for a reasonable doping fraction) the doping dependence. For instance, in [^J the max- 
imal length of the system studied was 16 and the doping was represented by just two holes. 
There are significant size effects and it is also difficult to extract directly information about 
correlation functions. DMRG approach allows to work with much longer ladders (up to 100 
cites in length), but to extract frequency dependent information is difficult. In any case, 
an analysis of numerical results frequently requires some input from field theory (see, for 
instance, (fl,!^,!^), and it would be highly desirable to have a description valid for unequal 
velocities. This is just what the available analytical techniques fail to provide. Although 
they have been very successful in describing the low energy sector of the model with equal 
velocities, it is not clear what happens if the velocities are different, as one may expect it to 
be in experimentally relevant situations. Another major difficulty is related to the richness of 
the fermionic ladder phase diagram which includes various phases with superconducting as 
well as Charge Density Wave quasi- long-range order (see [lOj and references therein). Since 
the analytic approaches start with unrealistic values of the bare interactions, parameters of 
the effective theory should be considered as phenomenological and one cannot, for example, 
trace their doping dependence. 

In this paper I present a procedure based on bosonization and subsequent refermionization 
of the 2-leg ladder Hamiltonian which allows one to derive an effective low energy field theory 
for a generic ratio of spin and charge velocities. My only assumption is that the resulting 
spectral gaps are small in comparison with the spin excitation bandwidth ~ J which seems 
to be consistent with the numerics Q]. The novelty of the approach lies in the judicious 
choice of the basis for the decoupled chains. Such choice greatly simplifies the form of the 
interaction and allows one to separate high and low energy degrees of freedom. Whenever 
is possible I compare my results with numerical ones. 

The paper is organized as follows. In Section II I discuss the model of doped ladder and 
derive the corresponding low energy effective field theory. In Section III I discuss different 
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regimes ( I call them underdoped, optimally doped and overdoped) with an emphasis on the 
underdoped one. In Section IV I discuss the problem of coupled ladders. There are separate 
sections for Conclusions and Acknowledgements. The paper has several appendices. 



II. FROM THE LATTICE MODEL TO THE LOW ENERGY FIELD THEORY 



As a starting point for our consideration I take the following model of coupled chains: 

+ V X2 (1) 



H = H X + H 2 - ^ t± [ C nl,a C n2,a + h.C 
n 

Hi = Yl { - 1 + h.c] + UN lt N tl } + IJ2 V < 

n 

V\2 = ^ |^Kim,± N ln N 2m + J nm Si ri S2 m 



N in Nir. 



nm 1 > m J ' im 



(2) 
(3) 



where N = ^2 a c^c a , S a = ^c+a^Ca' . In what follows I will assume that the interac- 
tions on chains are pretty much arbitrary and predominantly repulsive and the interactions 
between the chains (including the interchain tunneling t±) are small in comparison to the 
characteristic single chain energy scales. Below I will elaborate on these restrictions. The 
Hamiltonian ([T]) has U(l) xSU(2) xZ 2 symmetry. 



A. Derivation of the low energy field theory 



Away from half filling the low energy dynamics of the single chain problem is universal 
and is described by a sum of two Tomonaga-Luttinger models describing the charge and the 
spin sector. All details of the original lattice Hamiltonian ([2]) are encoded in few parameters 
such as the Luttinger parameters K C ,K S and spin and charge velocities v s ,v c . The SU(2) 
symmetry of the spin sector fixes the value of K s to be one. In the future we set K s = 1 and 
drop the subscript for K c so that K c = K. The Lagrangian of two noninteracting chains is 



r. 
r. 



dx 



(£ lc + C 2c ) + {C ls + C 2s ) 



2K 
1 



(4) 
(5) 
(6) 



v^\d T <$> ic ) 2 + v c {d x <S> ic ) 2 \, i = 1,2 
v;\d T $ is ) 2 + v c {d x <$> is ) 2 ] - 2nv s gJ R J L , 

z L 

where $j C , $ is are bosonic fields from the charge (c) and spin (s sectors) of z-th chain. The 
last term in fl6]) represents a marginal interaction of right- and left-moving spin currents JJj. L . 



These currents satisfy the SUi (2) Kac-Moody algebra; their explicit expressions in terms of 
bosonic fields are not necessary for the present discussion. For systems with predominantly 
repulsive interactions g > 0, g ~ 1 the current- current interaction is marginally irrelevant. 
The present universal description is valid below certain doping dependent cut-off A(<5) (5 is 
doping). For the Hubbard model with strong on site repulsion U » t the cut-off is 

A = min [ J « 4t 2 /U, At sin 2 (vr5/2)] . (7) 

The ratio v s /v c is also doping dependent since v c as a function of doping has a maximum 

-,,,,,,,,»«.,,„„,«-,,,-„„„, «,„-,.. — . 

My goal now is to choose a convenient basis of fields to treat the coupled chains. The 
first step is standard. Namely, I introduce bosonic fields 



$(±) = $(i) ± $(2) 



(8) 



together with their chiral components 0, = ($ ± 0)/2, B being the field dual to $. Let us 
consider, for instance, charge Lagrangian After transformation it becomes 

Cic + = ^ £ [^(dM^) 2 + v c {d x ^f] . (9) 

a=± 

Away from half filling the Umklapp processes are suppressed and field $c + ' ) decouples from 
the interactions (this statement is also supported by the direct calculations presented below). 
As far as the (c, — ) Lagrangian is concerned, I refermionize it: 

r + (d T - w c d x )r + l + (d T + w c d x )l + g K v c r + rl + l, (10) 

where 



r = — ^ e ^i~\ / = ^^e" i ^" ) , « 2 = 1, (11) 
V27rao V^ a o 

with k and ao being a coordinate independent Klein factor and a small distance cut-off 
respectively. The coupling gx is related to the Luttinger parameter K; for \K — 1| << 1 we 
have: 

g K ^ 2n(l/ K-l) (12) 



In view of the further developments I introduce Majorana fermions £ R L and tj Rj l related to 
the real and imaginary parts of conventional fermions r, I: 

r = {& + iVR)/V2, l = (e L + m)/V2, e + = e,V + = V- (13) 
In the bosonic language these Majorana operators are 

Vr = -p= sin(v / 4^0i" ) ), £| = -}L= cos(v / 4^0^ ) ). (14) 
y/irao y/Tra 

Majorana fermion fields can be expanded in terms of conventional creation and annihilation 
operators. For instance, we have 

V(x) = J2^ e ~ ikX + % +e H tik, <} = Sk P - (15) 

Substituting ( Fl3l) into ( fTOl) I obtain the Majorana fermion form of the Lagrangian for the 
antisymmetric charge mode: 

L(c-) = ^Vr(9t ~ v c d x )rj R + ^rj L (d T + v c d x )rj L + -£ 3 R (d T - v c d x )£ 3 R + 

~&(d T + v c d x )e L + 2g K v cVRVL e R e L - m 

Similar transformations in the spin sector (see [l^] for details) yield the following La- 
grangian density: 

1 3 

L » = 2 £ [*R& - ™> d M + X a L (d T + iv s d x ) X a L ] + Vex (17) 

a=0 

= -^^^(xki)(xki), (is) 

where Xi?L are a l so Majorana fermions defined as 

X°r = sin(v / 470i- ) ), xl = ~^== cos(V4^0«), 

= sin(v / 4^0W), x% = -£±= cos(v / 4^0W). (19) 

with k Sj -i- being new Klein factors. In the process of derivation I used the fact that 

4° + Jl a = ^OcXrXr, 4' a - Jl a = KrXr, (20) 



with similar expressions for the left currents. 



Now let us consider the interchain tunneling: 

Vt = t± ^2(ct ntU C 2n ,a + k.C.) « t± / Ax{^ R ^ 2 R,a + ^1£ 1<T 02£,<T + k.C.) + (21) 
n J 

where 

c n = e- ifc ^^(x) + e ikFn ij L (x), x = na (22) 

and the dots stand for interaction terms generated by virtual high energy processes. The 
latter terms contribute to the renormalization of interachain exchange and density-density 
interactions. Using for ipR^L bosonization formulae similar to (TTTj) I obtain: 

^ R ^2R,a + h.c=^T le^^e^^ - h.c] , (23) 

cr 

where Z is a nonuniversal amplitude which magnitude depends on the intra-chain interaction. 
Taking into account that the combinations of Klein factors k^k^ and commute with 

each other and their squares are equal to -1, we can choose them as equal to ±i and get for 

<m 

V t ~ 2it'j_ J dx T]rXr + VlX°l + interchain exchange, (24) 

where t' ± = Z 2 t±. In what follows I'll drop the prime superscript at t' ± . Recall that xo is 
one of the Majorana fermions from the spin sector. However, since it transforms as a singlet 
under action of the SU(2) group, the entire tunneling term is an SU(2) singlet, as it must 
be. Naturally, the tunneling entangles the charge and the spin sectors, but only the singlet 
fermion participates. As the next step I diagonalize the quadratic part of the Largangian 
( Tl6l l24l [ITl) containing r\ and x ■ Switching to the Hamiltonian formalism and dropping the 
superscript and subscript for x° R I get the following Hamiltonian for the right moving sector: 

H ' = Y1 [ v c kr lkVk + v s kxixk + t'±(xtVk + VkXk) , (25) 

fc>0 

with r^" 4 ", Xki Xk being conventional creation and annihilation operators. The spectrum 
is 



E±{k) = (v c + v s )k/2±J(v c -v s ) 2 k 2 /4 + t' ± z , k>0. (26) 
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Linearizing the spectrum close to the point Q = t' ± /y/v s v c where E_{Q) = and ignoring 
the upper branch of the spectrum which lays at high energies: E + > t' ± , I get 

2V C V S 



\k\«Q 



U 



v c + v s 



Vc 



V s + V c 
Vs ' 
V s + V c . 



1/2 



1/2 



2Vc lV2 



V s + V c 



(27) 



cos Qx^ R + sin Qx£ 2 R ) , (28) 



-lQx 



2v K 1 V2 



V s + V c 



- sinQx^ R + cosQx^|)(29) 



where 7?., 7?. + is a conventional (Dirac) right-moving fermion defined on the entire k axis 
and £( 1,2 ) are its Majorana components (its real and imaginary parts). For the left-moving 
fermions I get the same expressions with u —> —u and Q — > —Q. Substituting (129]) into ( IT6|) 
for energies less than t' ± I obtain the following Hamiltonian density: 



H c - + H s = H n + Vi 



inti 



(30) 



a=l,2 



0=1 



where Vi n t contains four-fermion terms. At this point it is worth to discuss briefly the 
problem of Fermi points. Naturally, one expects that the interchain tunneling splits the 
original Fermi point in two. From the above discussion one may get an impression that 
there is just one Fermi point. This is certainly not the case since even from Eqs. (I28I29|) one 
can see that there are oscillations at wave vectors kp ± Q. To check that no low energy 
modes got missing on the way, I compare the central charge of (|30|) with the central charge 
of theory of four noninteracting fermions. The latter one is equal to 4; central charge C=l 
goes into the symmetric charge mode (c, +) which decouples from the rest of the system. 
What remains is C=3 which coincides with the central charge of model of six Majorana 
fermions f )30|) . Indeed, every Majorana mode representing a half of conventional fermion 
carries central charge 1/2 so that 1/2x6 = 3. 

To obtain the interacting part of the Hamiltonian we have to take into account the four- 

Dosonized and refermionized 



fermion interaction in ( flOl) . as well as ( fl8l) and the properly 
interchain interaction. The latter one has the following form ( 



2l|) 



— cos(v / 4T$(- ) ) \(-J + V)J2 XrXI + (3 J + V) Xr X°l\ + V{d x ^f 
Tra L ^ 



o=l 



(31) 



Here J and V are related to the interchain exchange and density-density interaction and V 
to the density- density forward scattering. Since the forward scattering usually yields only 
weak corrections to coupling constants, I will put V = 0. In terms of the Majoranas we have 

1 



xlxl 



cos(v / 4^$^ ) ) = '-(vrvl + elei) = ^{ 



& & + If 



V c + V. 



As it was mentioned above, the symmetric charge mode $c + ' ) does not appear in the inter- 
actions and thus decouples. 

Bare fermions 

C+(l) C f (1) C f +(2) C f (2) C| + (l) C|(l) q + (2) C|(2) 



(32) 



Standard 



Bosonization 




SU(2) triplet 



U(1) doublet 



FIG. 1: The succession of transformations leading from the original formulation (1) to Eq. ([30 
The framed fields is the ones taking part in the final effective theory. 



Combining ([XBJ and ([32]) I obtain 



v mt = -2(e R e L ) [g.AxRxi) + 9c, ss mi - e R e L )\ - g P Aaki - e R e L ) 2 

3 3 

-2fc t (^-eM)E^)- 2 ^ E (x a RXl)(XRX b L ) (33) 



a=l 



a>b,a,b=l 
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where 
9c- 

9c,st 



{V-J)/2, g c , ss = 
(V -J) /2 + irgv c 



■ g K v s + V + 3J^)/2(v c + v s ), g p _ 



{v s v c ){V + 3J) 



[v c + v s 



9a,- 



irgv s 



(34) 



\v c + v s 

The complete effective Hamiltonian is a combination of the Luttinger Hamiltonian for $Jf 
and (13011331) . It has a similar form to the one obtained by the weak coupling bosonization 
procedure (see, for instance, lOj which notations I follow). There is one important difference: 
the velocities of different Majorana modes are now different. The resulting Hamiltonian, as 
expected, has the same Z 2 xU(l) xSU(2) symmetry as the original lattice model (1). In that 
sense its general form is dictated by the symmetry considerations. However, the microscopic 
derivation is not without use since it puts constraints on the model parameters which I am 
going to exploit. 



B. Some important operators. 

n 

Since the resulting field theory is qualitatively the same as the one described in [10J, all 
expressions for possible order parameters (OP) are the same as in 10] . Since I am going to 
concentrate on conditions for superconducting pairing, only two OPs are needed, one of them 
being the superconducting OP antisymmetric under chain permutation. In one dimension 
this property imitates (i-wave symmetry. Since a model of a single Majorana fermion is 
equivalent to the Quantum Ising model (see Appendix A, a more detailed discussion of this 
subject is given, for instance, in j^,[l5j|), it is convenient to express the order parameters 
in terms of order and disorder operators of the corresponding Ising models: a, /i (for x 
fermions) and S, M (for £ ones). So we have 



Asa ~ ^^{tTitTzaaSgSiMa - iM 3 Mx £2/1^3 }. (35) 

The other important operator is the 4fcp component of the particle density. Its impor- 
tance is related to the fact that the 4k p instability corresponding to crystallization of pairs 
competes with the <i-wave superconductivity. Meanwhile, the form of this operator remains 



controversial. 



. in 
H). 



10( it was conjectured on symmetry grounds, as in other publications (see, 



for instance, [8]). Here I provide a rigorous derivation. For the Hubbard model we have 

p(4k F , x) — F cos(4k F x + V8ir$ c ) (36) 
10 



where F is a non-universal amplitude. For the fermionic ladder one can define symmetric 
and antisymmetric combinations of charge densities: 



P 



(+) 



(4k F ,x) = 2Fcos(4k F x + VAti^ ( c +) ) cos(v / 4tt$^ ) 



iFcos(4^x + v^^){e R e L + T-xrrKki - e&ii} (37) 

and 

p {) {4k F ,x) = 2Fsm(4k F x + V4jr$ { c +) ) su^a/Itt^) = 
2Fsm(4k F x + y/4^$W)[riR& + fe] « 

iF^ s /(^ c + ^)sin(4fc F x + V47$( +) ) [e" i0a; (^2 + £|£ + ) - /i.e.] . (38) 

The latter formula indicates that weaker harmonics appear on wave vectors 4k F ±Q. 

The standard approach to fermionic ladders maintains that when the system scales to 
strong coupling and gaps emerge in various sectors, amplitudes of the corresponding OPs 
freeze and only the phase factors containing gapless charge mode $ c + ^ or its dual 0c remain. 
As a result the scaling dimensions of the OPs become much smaller than for noninteracting 
fermions leading to a strong enhancement of the corresponding susceptibilities. For instance, 
the SCd OP fl35|) being a bilinear combination of fermions; in the absence of interactions its 
scaling dimension d bare = 1, for decoupled chains dd eC ou P ied = 3/4 + 1/AK > 1. On the other 
hand, in the SCd phase of 2-leg ladder when the amplitude 

A = {a-irrarraSaSiMa - iMgM^^/is}. (39) 

acquires a nonzero ground state average the scaling dimension becomes 

dscd = 1/iK < 1. (40) 

One may think that 2k F Charge Density Wave OPs containing e^V^^i +) exponents with scal- 
ing dimension K/4 are more relevant when K < 1 (repulsive interactions). This, however, 
is not the case because it turns out that in the phase with nonzero SCd amplitude (1391) the 
amplitudes of all other OPs (except of the 4k F density one!) do not have nonzero ground 
state averages. The amplitudes of (|38|) also can never form because of the symmetry 
and ^ h L with a never pair). 
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III. DIFFERENT REGIMES. 



Although model (I30|33l) is not integrable, some general statements can be made concern- 
ing its spectrum, especially in the regions where two of the velocities v c and v s are very 
different. Since interactions (13"3"j) are marginal, for general value of parameters they can 
all flow either to weak or strong coupling. The numerics done on two-leg ladders usually 
joints to the strong coupling regime. However, a possibility of the weak coupling also exists 



271 ]. 28] and, as I will discuss further, apparently is realized in this model for v c ~ v s when 



one can apply the standard RG. 

As I am going to argue, the strong coupling regime is realized when the velocities are very 
different so that one excitation branch lays well above the other in most of the phase space. In 
that case I use the adiabatic approximation. Such approximation has been used extensively 
to study soliton excitations in the Peierls-Frohlich model pertaining to the polyacetylene 
problem where the slow subsystem is re pre sented by the lattice phonons. The representative 
capers on the model in question are 29j, the summary can be found in review articles 

mm- 



A. Underdoped regime v c « v s . Adiabatic approximation 



In the given case the fast particles are spin Majorana fermions Xa- I bosonize slow 
fermionic modes £*,£ 2 so that the Hamiltonian (I30|I33 j) becomes 



c= l - [u-\d T ef + u{d x ef] + x - [e R (d T - w c d x )e R + e L (d T + w c d x )e L 



A 



9s 



+ 



a=l 



iv s d x )xR + xli^r + iv s d x )xL + (4-fermion interaction) 



(41) 



9s,st 



cos[v47T0], 7 = Ag cM /g p _. 



Here it is assumed that _ > 0. In the limit u « v s field 9 becomes static [32|. The £ 3 
fermion being nonclassical object on this stage should be neglected. I will also neglect for 
time being the interaction between x-fermions. This interaction will renormalize the mass 
gap, but will not produce any qualitative changes. Then t 
Peierls-Frolich one discussed in the classic soliton theory 



re problem is equivalent to the 



29l,|30|,(3jJ. Namely, the field 
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theory problem is reduced to the solution of differential equations for the fast Majorana 
fermion wave functions supplemented by a self-consistency condition: 



Exe = 
A(z)/ 7 



W S T 3 — + A(x)t 2 

ax 



Xe, 



Xe t Xe 



(42) 
(43) 



where xe is a two component vector and f 2,3 are the Pauli matrices. It should be emphasized 
that though A-field is static, it is coordinate dependent. 

As was shown in 29[ the latter problem is exactly solvable. All excitations have spectral 
gaps. For the present purposes it will suffice to consider just a single soliton solution of 
(142 f 43 1) . The self-consistent solution of ( 142P3|) corresponds to A(x) = —k tanh[k (x — 
Xo)/v s ], where k$ depends on parameters of the model. Dirac equation (1421 with such 
potential has two kinds of solutions. One solution corresponds to E = and is localized at 
the soliton center: 



Xr 

Xi 



1 



1 



exp 



*A(y) 



dy 



'-Jo 





exp 



x Hv) 



dy 



.(44) 



fc >0 \ / \ / fc <0 

The total energy of this solution is finite and originates from the A 2 -term in (jUJ). Therefore 
this excitation has a spectral gap (I denote it Mb). Its quantum numbers are supplied by 
Majorana zero modes (jUJ) which also include £ 3 . In the limit u = the soliton does not 
move; at u ^ it acquires dispersion and in that sense it is a slow particle. Further down I 
will discuss these excitations in more detail. 

Another solution is a scattering state of a (anti) soliton and a massive particle of mass ko'- 



Xr{x) 



Xi{x) 



2\[L 
2^L 



1 + 



- kv s + ik tanh[A; (x — x ) / v s ] 
-kv s + ik tanh[ko(x — x )/v s ]~ 



(45) 



(these wave functions are normalized) with energy 



E{k) = ^Jk 2 + (kv s 



(46) 



As I have said, the spectral gap k Q for this excitation is established via self-consistency 
condition 0430 . This gap is different from Mb, as well as quantum numbers of this massive 
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particle do not coincide with quantum numbers of the slow soliton since its wave function is 
not localized at x — xq. Namely, the fast particle carries quantum numbers of the Majorana 
fermion x a and hence it is S=l neutral spin exciton. In order to treat quantum corrections to 
its gap I now reinstate the four-fermion interaction and write down the effective Hamiltonian 
for the renormalized fast modes: 



H 

T 



fast 



(47) 



0=1 



a a>b,a,b=l 

Since k is related to the amplitude of A, it follows from ( I4T1) that it is related to the vacuum 
average: 

k„ = -i(V + ^ - JXCeki - &D) ~ (J ^'l';;" n(A/Mg> Afa (49) 

The fermion-fermion interaction in (|47p leads to a substantial renormalization of the bare 
mass. The latter is determined from the following equation: 

k 



1 + § ln(A/m s ) ' 



(50) 



As it follows from ( 1491 the spin gap m s becomes smaller with the increase of v c (that is with 
the increase of doping). 

Naturally, in the presence of gapless charge mode <S>[ +) the spin exciton is incoherent. 
Since the staggered spin density operator includes a bosonic exponent of this mode, the 



18]. The resulting com- 



exciton is emitted together with a cascade of gapless excitations 
posite object is frequently called " magnon- hole-pair bound state", though, strictly speaking 
this is not a bound state, but a composite. T his p ictur e ag rees with the numerical results 

adder obtained in l^ . 



for 5 = 1/8, 1/6 2-leg 
respectively. Fig. 1 in 



17j and [19| for 2 x 24 and 2 x 12 ladders 



16| and Fig. 4 in [19|] show a sharp peak centered around 7r(l — 8) 



which gradually merges into the continuum at larger wave vectors. 

Now I return to the slow modes. As it has been explained above, these are solitons 
and antisolitons. As was first suggested in 33(, they receive their quantum numbers from 
the zero modes of Majorana fermions (f44l) coupled to them (in the present context the 



reader can also benefit from the discussion in 



10|). Such zero modes (£ 3 mode are included) 



compose spinor representations of the 0(6) group. The latter ones are equivalent to the 
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representations of SU(4) group with Young tableaux consisting of one box (particles) and 
a column of three boxes (antiparticles). Hence the solitons carry quantum numbers of 
electrons (that is transverse momentum and spin) modulo total charge. In other words, the 
slow excitations are quasiparticles with electric charge stripped into the (quasi) condensate, 
in short, Bogolyubov quasiparticles. I denote their spectral gap as Mb- As I have mentioned 
above, I do not see any reason to have different gaps for different transverse momenta as in 



16j, 19] • Since model f l30|33|) is not Lorentz invariant, a general form the spectrum is not 



fixed. At small momenta we have 



9 



E B (p) = M B + ^ + 0(p% (51) 

where M* ^ M B due to the lack of Lorentz invariance. To get an idea of the relative size 
of the effective mass M* I notice that all kinetic energy of a soliton comes from the bare 
kinetic term in (jUJ). Thus in the leading order in u/v s the correction to the energy of a 
moving soliton 8(x — vt) is 

V ^h0 2 -^~ M ^^ (52) 

where / ~ v s /M B is the solitons's size. As expected, the kink turns out to be a heavy particle: 
M* /M B ~ l/uv s . This feature provides consistency to the adiabatic approximation. Since 
it becomes energetically disadvantageous for zero energy bound states of slow Majoranas £ 3 
to "sit" on kinks moving with velocity v > v c , the region of solitons's stability in momentum 
space is restricted by the condition v < v c or \p\ < M*v c < Mb/v s . 

The last remaining problem is whether the underdoped regime may support d-wave su- 
perconducting pairing. As was explained in Sec. IIB, such pairing requires a formation of 
the proper OP amplitude and this, in turn, requires a proper structure of the vacuum. As 
follows from fl33|) . the relation between vacua of £ and \ fermions depends on the sign of 
coupling constants g c>ss ~ ~g K v s + V + 3J, g c , st ~ (V + irgv c - J) and g a _ = (V - J)/2. 
When all these constants are negative the masses of Xa and £1,^3 fermions have the same 
sign opposite to the mass of £ 2 . As a consequence the average (ci (72(73) forms simultane- 
ously with (Exil^X^) ( or (/i 1 // 2 /ia)with {M 1 T, 2 M 3 )) and hence the entire SCd amplitude 
fl35|) freezes (see the discussion in [lOj). Since V, J > the conditions to be fulfilled are 

(g K v s -V)/3>J>V + ngv c . (53) 
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I conclude this long line of argument with the following comments. 

A qualitative picture of the underdoped regime emerging from these calculations is as 
follows. First and foremost, when condition (j53p for the interchain exchange is fulfilled, 
the phase has quasi-long-range superconducting order with <i-wave symmetry. This requires 
sufficient strength of the exchange interaction which agrees with the previous expectations. 
The doping dependence of the spin gap m s ( 1491 and 150]) originating from the presence of v c 
in ( 149]) is a new result. 

In arrays of ladders SCd order competes with what is variously called Wigner Crystal of 
Pairs or 4k p Charge Density Wave. The latter phase is probably realized in the telephone 
number compound [26]. There is a possibility of m s changing sign before v c and v s become 
comparable, that is within the validity of the adiabatic approximation (this corresponds to a 
violation of J > V + ngv c condition). This would lead to a transition to Ikp Charge Density 
Wave phase |lo| . 

Third comment. Usually doping is associated with incommensurability. However, model 
(jUJ) looks like commensurate Peierls-Frolich model; the reason for that is that the only charge 
participating in the game is the relative one, and the way it couples to the spins reveals its Z2 
nature: it is just cos \/4ti9. Although the entire system is doped, and correlations are indeed 
incommensurate, all this is accounted for by the total charge field $c + ' ) which is decoupled 
from the rest of the system. The spin SU(2) symmetry is, of course, unbroken, so the spin 
Majorana mass bilinears can appear only as a scalar {xrXl)i an d a scalar can only couple 
to another scalar. The nontrivial conclusion is that as a result of the above the situation 
for a doped fermionic ladder is analogous to the commensurate Peierls-Frolich model with a 
real order parameter A. 

Forth comment. By using the results taken from the theory of Peierls-Frolich model I 
was able to demonstrate that at v c << v s the fermion ladder has at least three types of 
excitations. They have all been found in the numerical papers [l(J,[l7j], jl9]. There are 
"fast" (v s ~ t 2 /U) spin S=l magnetic excitons, gapless charge mode and heavy (slow) 
excitations. The spin excitations are triplet ones (spin S=l), that is they are qualitatively 
the same as the excitations of the undoped spin ladder, described by the Hamiltonian similar 
to (14T|) [14]. They have a relativistic-like spectrum: 

E s (k) = jm\ + (v s k) 2 . (54) 
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As I have said, single S=l is always emitted in conjunction with the charge gapless modes. 
The slow excitations also have a spectral gap, but their spectrum is not relativistic (15 ip. By 
their quantum numbers I identify them as Bogolyubov quasiparticles. It is instructive to 
compare the described excitation spectrum with the spectrum for the case v c ~ v s studied 
in {9 1, 11]. In the C1S0 phase this spectrum coincides with the one for the 0(6) symmetric 
Gross-Neveu (GN) model. It does contain Bogolyubov quasiparticles as kinks and antikinks. 
In that sense the difference in velocities is not important, as it should be expected for 
topological excitations. The 0(6) GN contains 6- fold degenerate multiplet of vector particles. 
These ones are local with respect to our Majorana fermions £, \. I have found that at 
v c « v s this multiplet is reduced to 3-fold degenerate multiplet of spin excitons Xa- Survival 
of the other three components including the famous Cooperons remains an open question. 

B. Single electron Green's functions for single ladder 

Some features of the single particle Green's function whose spectral density is directly 
measurable by ARPES, can be calculated in the strong coupling regime even if the model 
f)30|33p is not integrable. I perform this calculation for the right moving fermions. The 
decoupling of the total charge mode ©i + ^ from the rest of the Hamiltonian translates into 
factorization of the fermion creation and annihilation operators. Namely, the creation oper- 
ator operator for right-moving electron on j-th chain is given by 

^t,M = ^ c+) [ Z t«J QX ± Z t^~ iQX ] > V = 1,2), (55) 

where operators Z^ ap create excitations in the spin and parity sectors, index p = ± corre- 
sponds to the transverse momenta and 7r. The correlation function of the bosonic exponents 
is known: 

(e i^ + )(r,x) e i^ + )(0,0) ) = (r + ix /„ c )-l/4 L2 + ( X / Vc f~\ ^ K+1 ' K ' 2 ^ m (56) 

The excitation with smallest energy produced by Z + is a quasiparticle. It can also 
produce combinations of quasiparticles with magnetic excitons, but these processes require 
higher energy. Therefore in order to calculate the spectral function close to the low energy 
threshold Mb we need to determine just one matrix element corresponding to the emission 
of a quasiparticle. For Lorentz invariant models this can be done rigorously following the 
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standard procedure explained in [34], |35|. Although model (I30|l33j) is not Lorentz invariant, 
for energies not very different from the gap Mb when only slow solitons are excited this not 
be that important. In that case I will just treat (joTj) as the first terms in the expansion of 



a/ M b + (vp) 2 with v ~ y/v c v s . 

Taking these arguments into account I proceed further and assign operators with Lorentz 
spins. Recall that operator Oh with Lorentz spin h under Lorentz rotation 

vt = (t' v) cosh 9 + x' smh. 9 , x = (t' v) smh 9 + x' cosh 9 (57) 

transforms as 

O h (x,t) = e M O h (x',t'). (58) 

The right (left)-moving electron is a spinor and hence carries Lorentz spin 1/2 (-1/2), the 
chiral bosonic exponent carries Lorentz spin l/8(-l/8) and hence the right (left)-moving Z- 
operator must carry Lorentz spin 3/8 (-3/8). This means that the matrix element between 
the vacuum and quasiparticle excitation with energy e = Mb cosh 9 and momentum p = 
v^ 1 Mb smh. 9 parametrized by rapidity 9 is 

(9\Z+ aj± (x)\0) = Ze ±iw/4 e- 3e/8+iMB{x/v)s ' mhd . (59) 

where Z is a normalization constant and v is of the order of v c ~ u. The phase factor e ±l7r//4 
is necessary to maintain the (i-wave symmetry of the order parameter (see the discussion 
below). Substituting this expression to the Lehmann expansion for the Green's function 
and using (I5T)|) . I arrive to the following leading asymptotics of the single electron Green's 
function: 

(^(r,x)^+ ffy (0,0))~ (60) 
(r + ix/v c ) VV + (s/tt c )g](*+V*-2)/i6 Kv^ixTv) K H MWT2 + {X/V)2 



with v = ±1. Green's function (160]) is highly incoherent and is qualitatively similar to the 
Green's function of the 1/2-filled Hubbard model with different charge and spin velocities 



calculated in 



34J. For the sake of simplicity I restrict the calculations for the case v 



In that case for K — 1/2 we get 



.4 



u,p + vQ + k F ,k± = ?r(l + v)/2 ~ (61) 
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M 



-5/16 



(cu — vp) 



B 



— 33/32,9/32,5/16: 



M B + ( V P) 2 ~ ^ 

Ml 



(62) 



L 2 - (vp) 2 - Ml 

where v = ±1. The spectral function vanishes for u 2 < (vp) 2 + M B . The above expression is 
valid below the lesser of two energies: the three-particle threshold 3M# or the threshold for 
emission of two S=l excitons 2m s + Mb- It is also interesting to note that Q = t±/ \Jv c v s is 
non-universal and can be quite large. Therefore in the presence of interactions the transverse 
tunneling may shift the spectral weight very substantially from the single chain Fermi points 
Uf ± tj_/vp. 




FIG. 2: The spectral function A(u, ±Q + hp) = A(—u, ±Q + hp) (fBTj) as a function of 
uj/Mb > (blue). The sharp maximum at w = e(k) is replaced by a power law singularity. 
The magenta graph shows the same function convoluted with the Gaussian resolution function 
(3/^M B )exp[-9uj 2 /M 2 }. 

Two comments are in order. First, in the described energy range the spectral function 



is particle- hole symmetric which agrees with the early findings 18| and disagrees with the 



arguments by 23|. Particle- hole asymmetry is expected to appear at energies of order of 
the ultraviolet cut-off and is related to deviations of the bare band dispersion from linearity. 
Second, although the positions of the peaks in the spectral function depend on the transverse 
momenta, I have not found any difference in the spectrum for p y = 0, 7r which is reflected in 
the identical behavior of A(kp + Q) and A(kp — Q). This contradicts the numerical results 
of [l6]. 
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Although the anomalous Green's function vanishes for an infinite ladder, for finite ladder 
it may be nonzero provided proper boundary conditions are in place. This idea was exploited 



in 



19] 



20J to extract information about pairing. According to (j55p the anomalous Green's 



function is 



Fij(T,x) = ((lpR^,i(T, X^lpL^jiO, X 2 ))) = f(L,T,X 1 ,X 2 )J : ij{T,X 1 2) 

f(L,r,x 1 ,x 2 ) = ((e-iv^+)M) e iv^ +) (o,o)^ (63) 
T = (([z Rit+ (r,x)Z L>u+ (0,0) + [25^ - l]^ t ,_(r,x)^_(0,0)})). (64) 



As I said, the first factor in the last line of (154]) vanishes in the infinite system, but is finite 
in a finite one. It does not contain information about pairing and I do not discuss it. The 
information about pairing is contained in the second factor which leading asymptotics is 
determined by the formactors (1591) : 

1 



T~K (M b ^t 2 + (x/v) 2 ) 



[u + iO) 2 - (vp) 2 - M 



2 ' 

B 



(65) 



This expression is valid at frequencies less than 3Mg. In this region it looks like the standard 
Gorkov pairing function. The extensive incoherent part, found in [1_9|, 20] apparently appears 
above 3M# when the simple approach described here does not work. 



C. Optimal doping v c ~ v s 

n 

As it was established already in [9(, at v c ~ v s the interactions scale together either to 
strong coupling or to zero. In the former case the symmetry of the model dynamically in- 
creases at low energies and becomes U(l)xO(6). This regime has been thoroughly described 
in js], 11], 12]. It has been generally assumed that this situation is realized in t — J lad- 
ders. However, as it was pointed out in j^jj^sl, this is not necessarily the case. Indeed, if 
one takes as bare couplings ( 1341) and assumes that J > V (recall that this was one of the 
conditions for SCd in the underdoped regime) and vgx » V, J g ~ 1, the system will 
scale to the weak coupling phase C2S2. 
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Antinodal gap 



Critical point 



Doping, (-fc/ v s ) 



FIG. 3: The schematic phase diagram of doped 2-leg ladder. 



D. Overdoped regime v c >> v s 

Now the fast fermion mode is £ 3 . For this mode we can repeat the arguments presented 
in Sec. IIIA. Namely, the fast mode acquires a spectral gap (I denote it m 3 ); its spectrum is 



Its symmetry is Z 2 . There are also slow gapful solitons whose quantum numbers are the same 
as for the underdoped regime. Their spectrum is not relativistic. Beyond these statements 
I can provide no reliable information. 

IV. LADDERS BECOME STRIPES. 

The burning question is how properties of a single ladder are going to survive when one 
assembles ladders together. In this setting doped ladders ("rivers of charge") are separated 
by undoped ones. Such scenario have been considered many times before in the context of 
stripe phases and here the present approach differs only in details. 

The salient feature of the underdoped 2-leg ladder is a separation of scales between 
different excitations. Those ones include the gapless charge mode, the fast spin excitons 
with spin gap m s and the slow Bogolyubov quasiparticles with gap Mb- This separation 
is likely to survive in the bulk system, at least in some areas of the Brillouin zone. Since 
charge bosonic exponents enter into expressions for all relevant operators (such as 2k F and 
Akp components of spin and charge densities as well as interladder tunneling operators), 




(66) 
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transverse propagation of these excitations will be overdamped by emission of gapless phase 
fluctuations. The coherence will emerge below the 3D superconducting transition when the 
charge mode condenses (however, see the discussion of the single particle Green's function 
below). It is reasonable to assume that the overdamped Bogolyubov quasiparticles of our 
model correspond to the single electron spectrum in the antinode regions (see a more detailed 
discussion below). On the other hand, the direct interladder tunneling will give rise to Fermi 



pockets in the nodal direction, as was found in 



3- 



The precise position of the Fermi pockets 



at the Brillouin zone diagonals together with their orientation constitutes a difficulty for the 
present theory which it shares with all models of stripes (see, for instance, 

The spin excitons are probably responsible for the famous "hourglass" spectrum observed 



in the neutron scattering experiments 
different from the one performed in 



391 ] . The calculation presented below is not that 
where the doped ladder was described by the 0(6) 
Gross-Neveu model. Imagine doped 2-leg ladders running parallel with undoped ones so 
that a unit cell contains one doped and one undoped ladder. The most singular operator in 
the spin sector is the 2kp component of magnetization 

S_ = n cos(v / ^$^ _) ) cos(2A; F x + \/?r$i +) ) + n+ sin( v / ^$^ ) ) sin(2A; F x + V^ +) ) (67) 

Since the sector of $c is effectively frozen, we can consider all amplitudes to be static: 

S_ = Dn cos(2A; F x + y/n^), D = (coa{y/n^)). (68) 

So the spin operators in the underdoped regime are almost like the ones for undoped spin 
ladder, but slightly "softened" by the gapless charge mode. Then the effective interaction 
between the staggered magnetization of the undoped chains N and the magnetic modes of 
doped ladders is 

gN k (n_ k+7rS + n_fc-,rf) , (69) 

and the RPA spectrum is determined by the equation 

uj 2 -A 2 - (Vk) 2 - g 2 cos 2 q \— — L- - — + 1 — 1 = 0, (70) 

where q is wave vector transverse to the ladders and A is the spectral gap of the undoped 
ladders. When m s < A this equation yields the hourglass spectrum. 
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A. Superconducting fluctuations 



At energies smaller than the quasiparticle gap Mb and spin exciton gap m s one is left 
with the gapless charge excitations. In the phase where SCd order parameter acquires a 
finite amplitude, the 4k F component of density (157|) also acquires a finite amplitude. These 
operators from different chains couple together via Josephson coupling and the Coulomb 
interaction. This results in the following effective low energy Hamiltonian for the charge 
mode (I drop the superscript (+) and subscript c): 

Hcharge = ^ / dxU r + U C oul (71) 
r J 

Ucoui = 7r:2^ dxdx I, A 2 h m (72) 

271 ^ J y/(r-r') 2 + (x-x') 2 

K/(r, r') cos [V^(© r - ©r') - — / dlAl, (73) 



c 

where A is a vector potential, V c originates from 4k p component of the Coulomb interaction 
matrix element and V c is generated by the interladder Josephson coupling. The long range 
Coulomb interaction (172"]) is often omitted from consideration though it is important to 
get the right electrodynamics. Model ( ]7T]) is directly applicable to the telephone number 
compound. 

It is also instructive to write model (I7T|) in the Lagrangian form: 

N 

L = J2 dx £n, (74) 

n=l J 

£n = V -^-{d x § n ) 2 + id T e n d x $ n + ^{d x Q n f + U int , (75) 

where interaction term Ui nt contains all cosine terms. At K = 1/2 the two cosine terms 
in the interaction have the same scaling dimension and the competition between the super- 
conductivity and Wigner crystal ordering of pairs becomes strong. Recall that K = 1/2 
corresponds to the U » t limit of the Hubbard model and is a very realistic value in the 
present context. Then at Vj > V c the ground state is a superconductor and at Vj = V c there 



is a first order phase transition to the Wigner crystal state 40|. The latter state may emerge 
even at V c < Vj if the hopping of pairs between the ladders is suppressed by a magnetic 
field. Since field $ couples to disorder, In the presence of disorder the pairs localize. It is 
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well known that the suppression of superconductivity in underdoped cuprates reveals the 
state with weakly insulating properties 4l|, 42]. 

In Appendix B I consider exactly solvable model of two coupled ladders. 



B. Single particle Green's function in the striped phase 

In Section IIIB I studied the Green's function for a single ladder. Now I am going to 
generalize this calculation for a striped phase. Direct interladder tunneling produces bound 
states with quantum numbers of electrons. So, one may think that when ladders are coupled 
together the quasiparticles are reconstituted. However, since the tunneling matrix element 
is momentum dependent and changes sign throughout the Brillouin zone, the bound states 
(quasiparticles) are not created at all momenta, but only in those parts of the Brillouin zone 
where the tunneling process pushes the states below the continuum. It is well known that 
in the cuprates there are discernible quasiparticle excitations near the nodes and completely 



incoherent ones at the antinodes 



43j,|44|. Theoretically emergence of quasiparticle pockets 



in systems of weakly coupled chains and ladders have been discussed in |37|, [45] and more 



recently in 



46]. 



Below I ignore the direct tunneling which may be a good description for the antinodal 
region where the direct tunneling between the ladders does not create quasiparticles. As far 
as virtual tunneling processes are concerned, they generate the Josephson coupling between 
the ladders and are taken into account in effective Hamiltonian (17T1) . As for a single ladder 
case single electron excitation is a Bogolyubov quasiparticle dressed by fluctuations of the 
gapless charge mode Gc + ' ) . As for a single ladder one can assume that these two types of 
excitations are decoupled and the electron creation and annihilation operators are factor- 
ized as in (IB3|) . However, the charge mode is no longer one-dimensional which makes the 
calculation slightly nontrivial. 

For the sake of simplicity I will neglect the long distant Coulomb interaction fl72|) . In the 
cuprates this may be actually justified in the normal state where the Coulomb interaction 
is screened by nodal quasiparticles. Then in RPA the charge mode has the phonon-like 
spectrum 

n 2 = (v c q x ) 2 + (v/a) Ji sm 2 (qe t /2), (76) 
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where are the elementary lattice vectors. Above the transition there is a temperature 
dependent gap, but at the moment I ignore all temperature effects. In the absence of the 
long range Coulomb interaction model (I73p is Lorentz invariant in (r, x) space. Therefore 
one can still assign Lorentz spin to the chiral bosonic exponent 



h = e 1 ^^ (r) ( 77 ) 

and its value is still 1/8. Therefore the matrix element for emission of one bosonic excitation 
with transverse momentum and rapidity 9 defined as 

Q(9, qj_) = m(q±) cosh 9, q x v c = m(q±) sinh 9, m 2 = (v /a) Jj sin 2 (qe^/2) (78) 

i 

is given by 

(qx^le^^^lO) = ^(qje^e*/ 8 , (79) 

where Z is an unknown coefficient. Since in more than one dimension contribution of 
multiple emission processes is usually small due to the phase space factors, I will not consider 
such processes here. The contribution of the single emission to the correlation function of 
two the chiral exponents positioned on the same ladder is 

/ d D q ± Z( q± )( T ~ ^£ ) 1/8 K 1/4 [m(cu)^ + (x/v c )A (80) 

where D is the number of transverse dimensions. Simplifying matters again I take v c = v 
and get for the total correlation function 

G RR (T,x;r = r') = (81) 

, T u/v y/2 f d D Z ( )K \ m ( ^ T 2 + ( X / V y]K 3 jM B ^Tl + {X/Vf 

\t + lx/v/ J L J L J 

Assuming further that m ~ |gj_|, Z — const and D = 1 for u 2 > M B + (vp) 2 we get for the 
spectral function 

Au^,P + Q) ~ ^v^WFfn/S, 5/8, 1, *3 + %>'-"'' ) ( 82 ) 

It may seem counterintuitive that the increase of dimensionality of the phase fluctuations 
leads to a stronger suppression of singularity in the Green's function than for a single lad- 
der. The mechanism of this increase is related to the factorization of the electron operator 
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FIG. 4: The spectral function A(lu,Q) (l82j) as a function of oj/Mb (blue figure). The threshold 
singularity present for single ladder is removed by the emission of soft phase fluctuations (compare 
with Fig. ([2]))- The magenta figure represents the spectral function convoluted with the Gaussian 
resolution function exp[— 9u) 2 /M B ]. 



already discussed in Section IIIB. This factorization is similar in spiri ; 
factorization frequently used in gauge field theories of the cuprates (see 



; to the holon-spinon 
47] , |48] for a review) : 



4>a = hf a 



(83) 



where operator h carries charge and f a carries spin. The difference is that holons in that 
theories are pure bosons and hence can easily condense. The state with condensed holons is 
Fermi liquid where all incoherence is lost. There have been always qualitative arguments that 
the strong gauge field interaction between holons and spinons somehow prevents the former 
ones from condensation, but these claims have never been substantiated by calculations. In 
the present case there is no such danger. The role of holon is played by the bosonic exponent 
( I77|) which has a nonzero Lorentz spin and therefore cannot acquire nonzero average. On 
the other hand, an object with fractional Lorentz spin becomes very incoherent in more than 
one dimension which explains the effect. 
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V. CONCLUSIONS 



The doped 2- leg ladder demonstrates the most salient features of the bulk cuprate mate- 
rials. 

• When the interchain exchange dominates over the interchange Coulomb repulsion the 
phase with one-dimensional analogue of <i-wave superconducting order parameter is 
formed. It is understood that for a one-dimensional ladder this means only quasi long 
range order. 

• When doping increases the system goes through different regimes distinguished by their 
excitation spectra and behavior of correlation functions. They have many properties 
similar to underdoped, optimally doped and overdoped regimes of the bulk cuprates. 

• The excitation spectrum contains both Bogolyubov quasiparticles and collective 
modes. In the underdoped regime energy scales of all excitations are well separated. 
Collective modes include spin-1 neutral excitons and superconducting phase fluctua- 
tions. 

• All excitations are incoherent being overdamped by gapless phase fluctuations. Hence 
the single particle spectral function does not exhibit a sharp peak, but a broad maxi- 
mum (pseudogap phenomenon). 

• The ratio between spectral gaps of quasiparticles and spin excitons is interaction de- 
pendent. The spin exciton gap decreases with doping. 

• The magnitude of all gaps decreases with an increase of doping. There is a possibility 
that the optimally doped regime where charge and spin velocities are approximately 
equal is completely gapless. 

When ladders are assembled into a quasi-two-dimensional array 

• Single particle spectral weight becomes even more incoherent when ladders are assem- 
bled in an array and superconducting phase fluctuations cease to be one-dimensional 
(see the discussion in the end of the previous Section). The long incoherent tail of the 
spectral function (see Fig. 4) is its an inherent property and not a background. 
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• Likewise, a famous hourglass spectrum emerges in the dynamical spin susceptibility 
of the stripe phase. 

• Suppression of superconducting coherence by magnetic field drives the system into the 
insulating state with localized pairs. 

As I have discussed above some of these properties have been detected in the numerical 
calculations. 
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Appendix A: Bosonization rules for Ising operators 

The Quantum Ising model 

H = E - h <) ( A1 ) 

n 

is equivalent to the model of a single Majorana fermion. In the continuum limit its spectrum 
is relativistic e(k) = \/k 2 + m 2 where m — 1 — h. Since two Majoranas constitute one 
conventional fermion, two Ising models are equivalent to a model of a single Dirac fermion 
of mass m. The latter model can be bosonized and is equivalent to the sine Gordon model. 
Using this chain of equivalencies one can establish a correspondence between order and 
disorder parameters of the two Ising models and operators of the Gaussian model. The 
order parameter operator a is defined as a continuum limit of a*, its dual operator /i is a 
continuum limit of 

oo 
j=n 
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The bosonization rules for Ising model operators £1X6 cLS follows. 



cos(a/7T$) = fM\fJ,2, sin(v / vr$) = 
cos(\/vi : 0) = Mi°"2) sin(-v/7T0) = 



(A3) 



The above order parameters are sensitive to sign of m; (cr) 7^ for m > and (//) ^ for 



m < 0. This property is used in the main text of the paper. I refer the reader to 
further discussion of applications of Majorana fermions. 



15] for a 



Appendix B: Phase fluctuations in 4-leg ladder. Exact results. 

To illustrate a competition between Wigner crystallization of pairs and superconductivity, 
I consider a model of two coupled ladders. I define the modes 

$i )2 = ($+ ± $-)/v / 2 (Bl) 

and the symmetric mode (+) decouples. The Hamiltonian for the asymmetric mode is 



(d x $-Y + (d x Q. 



+ V c cos 



47T$_ 



Vj cos 



4tt0_ 



B2) 



where K = 2K and 



V c ~ U{4k F )M B , Vj ~ {t 2 /M B ){M B m t 



n3/8 



(B3) 



After refermionization I get the sum of two off-critical Ising models: 



iv r . 



H- = -T^i-PRdxPR + Phd x pL - Vr9xVr + VLd x 7] L ) + 
4nv c (K - 1)prPlVrVl + 2i(V c + Vj)p R p L + 2i(V c - Vj)rj R r] L . 



(B4) 



Consider K « 1/2. Neglecting the four-fermion interaction I obtain from (1B4I) that the 
amplitude of the superconducting order parameter diminishes when Vj approaches V c from 
above and vanishes at Vj < V c : 



giv^e _ e iv^e +e i0fe_ _ e i^e+ 



/iicr 2 + \o\\li 



(B5) 



The amplitude is 



2M/ 



(^a 2 ) ~ {Vj - V c 2 ) 



(B6) 
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At V c > Vj one gets the same expression for the amplitude of the 4kp density wave. 
Magnetic field adds to (IB4I) the term 



2ih(r] L p L - rj R p R ) 



(B7) 



The spectrum becomes 




{v c kf + h 2 + (vj + v 2 ) ± A{v c kh) 2 + 4h 2 vJ + v 2 v u 



(B8) 



The critical point where u_(k = 0) = is achieved at 



h 2 = (V 2 - V 2 )/A. 



(B9) 
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